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Abstract 

We calculate the baryon asymmetry produced at the electroweak phase transition by quasi- 
degenerate third generation sfermions in the minimal supersymmetric extension of the Standard 
Model. We evaluate constraints from Higgs searches, from collider searches for supersymmetric 
particles, and from null searches for the permanent electric dipole moment (EDM) of the electron, 
of the neutron and of atoms. We find that resonant sfermion sources can in principle provide a 
large enough baryon asymmetry in various corners of the sfermion parameter space, and we focus, 
in particular, on the case of large tan /3, where third-generation down-type (s)fermions become rel- 
evant. We show that in the case of stop and sbottom sources, the viable parameter space is ruled 
out by constraints from the non-observation of the Mercury EDM. We introduce a new class of CP 
violating sources, quasi-degenerate staus, that escapes current EDM constraints while providing 
large enough net chiral currents to achieve successful "slept on-mediated" electroweak baryogenesis. 
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I. INTRODUCTION 



The origin of the matter- antimatter asymmetry observed in the universe [haryogenesis) 
is one of the greatest puzzles of modern cosmology. No "standard model" of baryogenesis 
exists, nor does it seem possible to accommodate baryogenesis within the Standard Model 



of particle physics [1], although numerous hypothetical frameworks 



lave been investigated 
for reviews). 



within various extensions to the Standard Model (see e.g. Refs. 

m 

Of all the known options for a dynamical generation of the matter-antimatter asymme- 
try, one that may soon be conclusively tested is the framework of electroweak baryogenesis 
(EWB). This scenario generically entails the existence of particles that have electroweak 
quantum numbers and that are light enough to be in thermal equilibrium at the electroweak 
phase transition, when the universe had a temperature of around 100 GeV. A host of phe- 
nomenological consequences then follows: 

(i) light, electroweakly-interacting particles should be produced at the Large Hadron 
Collider (LHC), making EWB testable with colliders (see e.g. Ref. ^, 5]); 

(ii) large enough CP violation and light enough new electroweakly-interacting particles 
generically produce sizable electric dipole moments (EDMs), via one- or two-loop effects, 
to the level probed, or close to being probed, by experiments (for recent studies see e.g. 
Ref. 0,9); 

(iii) if the lightest electroweakly-interacting particle is stable, it could be the dark matter; 



recent studies j8l-ll0[| indicate that if this is the case, the resulting dark matter particle should 
be detectable with the current generation of dark matter search experiments; 

(iv) bubble collisions at the necessarily strongly first-order EW phase transition produce 
gravity waves with characteristic frequency /ew ^ T^w/Mp^ which, today, have a frequency 



/ew " ^ 0.01 mHz, in t 



peaks (see e.g. Ref. [U 



le range where the sensitivity of planned space-based interferometers 



Broadly, the above phenomenological probes make EWB a testable scenario, even though 
quantitative predictions necessarily depend on the specifics of the model under considera- 
tion. For several reasons, independent of baryogenesis, a supersymmetric extension to the 
Standard Model has long been advocated as a particularly compelling one, with the minimal 
supersymmetric extension (MSSM) its simplest paradigm. 

The MSSM contains enough degrees of freedom to overcome the two key difficulties 



encountered in embedding EWB within the Standard Modeh (i) producing a strongly first 
order phase transition to prevent baryon number wash-out, and (ii) providing large enough 
sources of CP violation to quantitatively explain the observed baryon asymmetry. The first 



difficulty is overcome either by advocating a light sea" 
numerous studies in this vein, from Ref. jl^ to Ref. js]). 



ar top (hereafter "stop") (there are 
by invoking a non-minimal Higgs 



sector comprising additional singlet fields (see e.g. Ref. 115, ll6| and references therein), or 
by advocating higher-dimensional operators (e.g. Ref. IITI] and references therein). With 
regards to (ii), additional CP violating sources are also easily found in the MSSM parameter 
space. 

Currently, the electroweak scale is under close scrutiny at the LHC, and EWB is either on 
the brink of discovery or of being convincingly ruled out, at least in its simplest incarnations. 
To address the question of whether or not EWB might be nature's mechanism of choice to 
produce the baryon asymmetry of the Universe (BAU), it is imperative to investigate in 
great detail all possible open corners of the MSSM parameter space that could be related 
to this scenario. Recent important progress has been made in evaluating the impact of the 



LHC program on point (i) above, with strong implications for the light-stop scenario jl8| but 
also for frameworks with non-minimal Higgs sectors jlO^. In the present study, we concern 
ourselves, instead, with issue (ii) above: which MSSM sector could have provided the CP 
violating sources necessary to feed the production of a baryon asymmetry at the electroweak 
phase transition? We leave aside here the (equally interesting, but separate) question of 
the mechanism responsible for producing a strongly enough first-order electroweak phase 
transition. The present study is intended to provide a qualitative and quantitative step 
forward in closing in on the allowed parameter space for supersymmetric EWB models. 
Supersymmetric EWB has been extensively studied in many different contexts in the 



recent and not-so-recent past (see e.g. Refs. [4, [sl, Isl, l20l-l33[|). Most previous studies of 
EWB in the MSSM have insisted and focused on the CP violation relevant to produce the 
baryon asymmetry as occurring in the higgsino-gaugino sector. This scenario generically 
forces a close degeneracy between one (or both) of the electroweak gaugino soft supersym- 
metry breaking masses Mi, M2, and the supersymmetric /x parameter, perhaps raising con- 
cerns about naturalness since Mi^2 and /x have entirely disconnected origins. Additionally, 
electroweak two-loop contributions to EDMs (see Ref. j3] for the complete set of relevant 
two-loop contributions) provide stringent constraints on the viable parameter space. 



Perhaps more importantly, there does not need to be any CP violation in the gaugino- 
higgsino "fermionic" supersymmetric sector for EWB to be successful. As noted long ago Q], 
CP violation in the stop sector is potentially an important source for EWB in the MSSM (and 
beyond IssI). Stop sources have not received much attention in previous studies primarily 
because of the superficial tension that typically exists between producing the correct BAU 
and having a phenomenology consistent with observation. This point can be illustrated in 
short as follows: the phase transition in the MSSM is typically made strongly first order by 
a light right-handed stop (typically in the sub- 100 GeV range - for an updated analysis see 
Ref. 15]); since the stop CP- violating sources are resonant for m~-^ ^ m^^, the parameter 
space most likely to produce the observed BAU has light, nearly degenerate stops. However, 
light stops are tightly constrained by the Higgs mass, elect roweak precision observables, 
and electric dipole moments (EDMs), along with null results from collider searches, thus 
making it difficult to have successful stop-sourced EWB along with a strongly first-order 
phase transition. 

If the ratio tan /3 of the two Higgs vacuum expectation values in the MSSM is sufficiently 
large, the Yukawa couplings of the down- type third generation fermions are generically of 
the same order as the top Yukawa coupling, therefore non-negligibly interacting with the 
Higgs fields. If this is the case, one would expect new contributions to the CP violating 
sources responsible for the generation of a net charge density of left-handed fermions, in turn 
:"ueling baryon number density production via weak sphalerons at the EW phase transition 



221]. Unfortunately, down- type fermion superfields do not contribute to the strength of 
the EW phase transition, since the cubic term they produce in the effective potential is 
counter-balanced by a l/tan/3 factor, making down-type contributions smaller than the top 
contributions by factors of m^^^/mt. Nevertheless, CP sources relevant for EWB from scalar 
bottom and r are generically significant in the large tan/3 regime. In this work, for the 
first time, we quantitatively study these sources, and we assess the impact of EDMs if CP 
violation is indeed sourced in the down- type (s) fermion sector. 

This manuscript is structured as follows: we first elucidate, in Sec. [Ill the physical model- 
ing of the electroweak phase transition relevant for the calculation of the baryon asymmetry 
of the universe (BAU), including details of the approximations we make and of the asso- 
ciated expected systematic effects. We then focus on stop sources in Sec. ITTTl and study 
several slices of the relevant parameter space. We consider constraints from Higgs searches. 




4 



EW precision observables, direct searches for supersymmetric particles with the LHC, and 
we devote special attention to constraints from EDM searches. Sec. [IV] then studies the 
large tan /3 regime and the new class of CP violating sources introduced in this study: quasi- 
degenerate sbottoms and staus. As for stops sources, we investigate all relevant constraints, 
with special emphasis on EDM searches. 

The main result of the present study is that the current EDM search limits - particularly 
the one obtained for the Mercury atom [3j] - eliminate the possibility that CP violating 
sources stemming from light and/or quasi-degenerate stops or sbottoms could be the main 
triggers for successful EWB. On the other hand, a new class of CP violating sources asso- 
ciated with third generation sleptons is subject to considerably weaker EDM constraints. 
Consequently, these sources can successfully generate the production of the net left-handed 
chiral charge needed to produce the observed baryon asymmetry at the EW phase transi- 
tion. This "stau-mediated EWB" is possible, however, only in relatively narrow strips of 
parameter space where the two stau mass eigenstates are almost degenerate. 

In view of the fast pace with which the LHC is exploring the electroweak scale, and es- 
pecially the supersymmetric sfermion sector, evaluating the potential relevance of sfermions 
to produce the observed baryon asymmetry appears to us as a timely topic. In addition, the 
program of searches for EDMs at the "intensity frontier" is also here demonstrated to be 
highly synergistic to the collider "energy frontier" ; searches for the EDM of multiple particle 
and atomic species is also crucial to testing the EW route to baryogenesis. Finally, while it 
is too early to draw strong conclusions, LHC searches for the Higgs might indirectly point 
to a scenario with light staus (see e.g. jsS]), potentially making the new source class we 
discovered all the more appealing. 



II. THE BARYON ASYMMETRY IN SUPERSYMMETRIC ELECTROWEAK 
BARYOGENESIS 

In the framework of supersymmetric EWB, the baryon asymmetry of the universe is 
produced by SU{2) sphalerons acting on a charge density of left-handed fermions, generated 
by CP- violating sources Sf^^ associated with the electroweak phase transition. As bubbles 
of broken electroweak symmetry nucleate and expand, the CP- violating phases between the 
supersymmetric particles and the background Higgs fields lead to the production of net 
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charge densities when (s)particles scatter off of the EWPT bubble waU. In the present 
study we concern ourselves with scalar sources associated with third-generation (s) quarks 
and (s)leptons, as their Yukawa couplings are much larger than those of their first- and 
second-generation counterparts. In addition, we focus on the large tan /3 regime, where third- 
generation down-type Yukawa couplings become comparable to the top Yukawa coupling, 
and therefore relevant in scattering oS of Higgs fields. The relevant part of the MSSM 
Lagrangian describing the associated CP- violating interactions in the gauge eigenstate basis 
reads: 

C D VttMAHl - ^i*Hf ) + ythbUAH', - ^i*Hl*) (1) 
+ VrfLfniArH', - ii*Hl*) - hHlH', + h.c, 

where CP- violating phases can arise between the various triscalar couplings Af^ fi^ and 
the Higgs soft mass parameter b. We henceforth denote this phase for species / as 0/ = 
Arg(/X74j6*). Without loss of generality, we will assume b to be real so that 0/ = ATg{fiAf) 
in what follows. 

In addition to the CP-violating sources, there are several CP-conserving processes arising 
from Eq. ([T]) that affect particle number rii for the relevant species in the MSSM. There are 
relaxation terms associated with chirality-changing particle scattering off of the Higgs vevs, 
with corresponding thermally- averaged rates Pf^. There are triscalar and supersymmetric 
Yukawa interactions given by Eq. ([T]) without replacing H^^ by their vevs; as discussed 
below, the assumption of supergauge equilibrium allows us to combine the rates for both 
types of processes which we write as P^^^. For the squarks, there are also SU{3) sphalerons, 
with rate P55, that produce 1st- and 2nd-generation squarks from a 3rd-generation density 
and vice- versa. Finally, weak sphalerons ultimately convert the left-handed particle density 

to a net baryon asymmetry with rate T^g- A complete set of expressions for these various 



sources can be found in Refs. 
the calculation. 
We follow Refs. 



23|, 1321], to which we refer the Reader for additional details of 



21 



25 



32| and work in the Higgs vev-insertion approximation, in which 
it is assumed that the sources in the bubble wall are strongest near the unbroken phase 
and where one uses a basis of SU{2) gauge eigenstates, expanding about fiavor- diagonal 
states in the bubble wall. This approximation tends to overestimate the resulting baryon 
asymmetry and clearly breaks down farther inside the wall where fiavor mixing cannot be 
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neglected. However the vev-insertion approximation is expected to characterize the produc- 
tion of the BAU to order unity accuracy [32]. Recent studies have worked out the flavor 



36| 



osciUations in the bubble wall beyond the vev-insertion approximation for a toy model 
and found qualitatively similar results to those obtained in the vev-insertion approximation, 
including a resonance in the various sources. Although a treatment beyond the vev-insertion 
approximation is desirable for an accurate assessment of EWB in the MSSM, since we will 
be looking at the baryon asymmetry across a wide range of parameter space with other 
inherent uncertainties, we content ourselves with the vev-insertion approximation, deferring 
a more detailed analysis including flavor-mixing effects to future study. Our results can thus 
be interpreted as a "best case scenario" for EWB with scalar sources in the MSSM, albeit 
we also show results that would correspond to a factor 10 smaller net BAU, to guide the 
reader to a more conservative interpretation. 

Proceeding within the outlined framework for computing the baryon-to-entropy ratio Yg, 
the weak sphaleron rate T^s is typically much slower than the rates for the creation and 
diffusion of the left-handed charge density Ul ahead of the EWPT bubble wall. This allows 
us to consider separately the diffusion equations for the various (s) particle densities and the 
creation of the baryon density n^, which is given, in terms of F^^, n^, and the bubble wall 
velocity as [20]: 

riB = / dz nL{z)e , (2) 

"^w J —oo 

where z is the distance from the bubble wall in the wall rest frame (neglecting the wall 
curvature) and where the unbroken EW phase corresponds to z < 0. The left handed charge 
density til given by the sum of the charge densities of the various left-handed quarks and 
leptons riL = + ^i) where the sum runs over all colors and generations and g^, k denote 
the difference of particle and antiparticle densities for each species. The charge densities 
entering into the expression for til are obtained from a set of coupled quantum Boltzmann 
equations (described below) which, once solved, allow one to compute Ub via Eq. ([2]). 

Detailed derivations of the quantum Boltzmann equations (QBEs) governing the gener- 
ation of the BAU have been discussed at length in the existing literature (see e.g. Ref. j32| 
for a full treatment) so we do not reproduce them here; in what follows we use the simplifled 
form of the QBEs discussed in Ref. [23], with some modiflcations. For each particle species 
we can deflne a corresponding chemical potential /i^, which is the fundamental quantity 



entering into the Boltzmann equations, related to its corresponding charge density by 

rii = ^kiiii + O (^^) , (3) 

where we have expanded in /x/T and the statistical weight for the density ki is given by 

6 r , xe^ 



h = g^— / dx Jx^ - mf/T^. (4) 

Additionally, as we will see in the following sections, for the cases we consider, the so-called 
supergauge rates, which drive chemical equilibrium between particles and their superpart- 
ners /x^ ^ /x^, are typically faster than the corresponding diffusion timescale Tdiff, defined 
in terms of the various diffusion constants and fc-factors in Ref. [32^. As a result of this 
"superequilibrium" condition, one can define common charge densities for the various par- 
ticles and their corresponding superpartners: Ui for right-handed up-type (s)quarks, Di for 
left-handed down- type (s) quarks, Qi for left-handed (s) quarks, H for the combined Higgs- 
Higgsino density, Ri for the right handed (s)leptons, and Li for left-handed (s)leptons (here 
i is a generational index). We also use the notation Q = Q^^ T = Us^ B = D^^ L = L^^ and 
R = R^^ while the fc-factors for these densities are defined by kj = ki + k~-. In terms of these 
definitions, the fermionic part of the density / (the quantity entering the weak sphaleron 
equation for the LH densities) is given by 

= (5) 
ki 

and the LH fermionic charge density is 

Two more observations allow us to reduce the number of equations needed to solve for the 
various densities. First, since weak sphalerons are decoupled from the Boltzmann equa- 
tions, baryon and lepton number are approximately locally conserved, so that the sum 
of all the densities vanishes at a given spacetime point. Second, since the first and sec- 
ond generation Yukawa couplings are negligible compared to corresponding couplings for 
the third generation, a first and second generation quark charge can arise only through 
strong sphalerons, and thus all corresponding charges will be produced in equal number, 
i.e. Qi = Q2 = — 2[/i —2U2 —2Di —2D2. Combined, these two relations imply 

8 



B = —{T + Q) 80 that we can eliminate the set of equations governing the B density as weU 
as aU of the other first and second generation (s)quark densities besides Qi. 

Given the above assumptions, the relevant set of Boltzmann equations to consider are: 
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(9) 



(10) 



(11) 



(12) 



We solve these equations in the so-called diffusion approximation, in which one introduces 
a diffusion constant for each species and assumes = DiVrii. The diffusion constants 
we use are those found in Ref. Dq = Dt = Dq^ ^ 6/T, Dh ^ lOO/T, Dl ^ lOO/T, 
Dr 380/r where T is the EWPT temperature, assumed to be 100 GeV. Note that the 
left- and right-handed (s)lepton diffusion constants are different; this is because of the SU{2) 
interactions active in the plasma for LH-densities. We neglect this difference for the (s)quark 
diffusion constants since Dq^t,q^ are determined primarily by SU (3) interactions which are 
non-chiral. 

With our framework in place, we can now compute the various sources and rates based 



on previous work in Refs. 



21 



25 



32| for the stop, sbottom, and stau cases. We assume the 



transition temperature Tc = 100 GeV throughout. 
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III. STOP SOURCES 



We begin our study by focusing on the scenario in which the observed baryon asymmetry 
originated primarily from the stop sector. When the stops scatter off of the spacetime- 
dependent Higgs vevs in the bubble wall, the CP-violating phase (pt arising between the 
tri-scalar coupling At and /x results in a non-zero expectation value of the current density 

in and in front of the wall, governed by 

d,i''^{x) = Si^{x,{n,}). (13) 

Here Si^{x^ {^i}) contains both the CP- violating source term as well as the CP-conserving 
chirality changing rates that also arise from stop scattering off of the Higgs vevs in the 
plasma, Yukawa interactions and strong sphaleron rates. To obtain the Boltzmann equa- 
tions for the stop case as in Eqs. fl7lfT2l). one must verify that the supergauge interactions 
governing the various particle and sparticle densities are in fact in equilibrium for the range 
of parameters we consider. Since we will vary the soft breaking masses of both stops, one 
should be concerned that in some regions of the parameter space, the supergauge interactions 
involving Ir and II will be slow (since these rates are Boltzmann suppressed) or kinemat- 
ically forbidden, since these rates arise from three-body interactions of the (s) quarks with 
gauginos. We plot the supergauge equilibration time scale in Fig. [1] along with rdiff for 
comparison (note that the kinematically forbidden region depends on the precise choice for 
the gaugino masses) . While the supergauge rate Vy is kinematically forbidden for very light 
RH stops, Ty is nowhere forbidden. This is because the latter is a sum of both W and B 
contributions which are disallowed for different Mg^ , and so when W interactions are disal- 
lowed, B interactions can still be active and vice versa (again, this depends on the details 
of the gaugino masses). Everywhere else both LH and RH rates are quite fast compared 
to Tdiff. The only other exception is the region corresponding to heavy squarks, where the 
baryon asymmetry is also expected to be suppressed. Since r^q <C Tdiff for most of the pa- 
rameter space relevant for EWB we work under the simplifying approximation that stop-top 
supergauge equilibrium holds in all regions of interest when computing Yg. Additionally, 
there are supergauge rates involving the other charge densities occurring in Eqs. fl7lfT2l): 
we have verified that the corresponding rates for Higgs and higgsino densities are also fast 
compared to rdiff for our choices of parameters, detailed below. The supergauge rates for 
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FIG. 1. Supergauge equilibration time scales for the RH (s)tops (Left) and LH (s)quarks (Right), 
where Mq^ (M^J = 1000 GeV in computing the RH (LH) stop rates and Mi = 100 GeV, Ms = 200 
GeV. Also shown is the diffusion time-scale Tdiff in both cases. The superequilibrium timescale is 
longer than Tdiff only in kinematically forbidden regions and for heavy squarks, where the baryon 
asymmetry is suppressed. 

the heavy squarks we consider are suppressed, and their equations decouple from the full 
set in Ref. [32]. As a result, the density Qi consists entirely of fermions, Qi = qi. 



In computing the baryon asymmetry we use the form of the sources computed in Ref. [21 1 
and related work, which exhibit resonant behavior for nearly degenerate RH- and LH- stop 
masses. We quote the form of the stop CP- violating source S-'^^ here, to allow straightfor- 
ward generalizations to the cases of sbottom and stau sources in the following sections: 



27r2 

X / Im 



f 

Jo 



nB{£*R) - riBiSi) _^ riBiSR) + nBiSi) 



(14) 



The various quantities involved are given by 

£l,r = ^L,R - i^L,R (15) 



uJl,r = ylkl +ml^^ (16) 
^^(^) = - (ex/T_i)2 (17) 
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where are the thermal widths of the stops which are (9(10~^T), are the spacetime- 
dependent Higgs vevs, v'^ = v^^ + v'^^ tan/3 = Vu/vd^ yt is the top Yukawa couphng, T is the 
EWPT temperature, Nc is the number of colors, and m^^ ^ are the effective stop masses given 
in terms of the corresponding soft breaking and thermal masses by ^ = ^ +M^.j^ ^. 
The dependence of Eq. (fT^ on the CP- violating phase (/)t is apparent. Also, the CP- violating 
source (and the chirality- changing CP-conserving rates) are manifestly spacetime-dependent, 
as they are proportional to the Higgs vevs. We use a simplified step-function profile for these 
rates and sources, deferring a careful treatment of the bubble profile to future study. Note 
also that we have omitted a temperature-independent contribution to the numerator of the 
second term in Eq. f fT4l) that appears in the corresponding expression in Ref. ^]. The cur- 
rent density from which the CP- violating source is derived must be properly normal-ordered 
through a subtraction of the zero-temperature matrix element. Implementing this normal 
ordering eflFectively removes the temperature-independent contribution to the numeratoj^. 

In addition to the CP- violating source Eq. f fMl ). we use the form of the relaxation, Yukawa, 
triscalar, and strong sphaleron rates worked out in Ref. j32|; we do not reproduce them here 
for brevity. Since we are interested in only the stop CP-violating source contribution to 
the BAU, we take the RH sbottom and RH, LH stau masses to be heavy which allows us 
to neglect the (s)bottom, (s)tau Yukawa rates - with heavy superpartners, only SM-like 
Yukawa interactions contribute to these rates, resulting in P^^^^r <^ ^yt^ ^ss in virtually all 
of the parameter space we consider. With this choice of spectrum we can also neglect the 
CP-conserving chirality changing rates P^^ for the (s)bottoms, which are suppressed by 
factors of (^^^ cot^/3 with respect to P^^ [23]. With these simplifications, the only source 
for B charge density are strong sphalerons, implying that B = —{Q + T) = —\Qi and 
consequently simplifying Eq. (|6j) to 

nL = ^g + 4(g + T). (19) 

Due to the relation between B and Qi and the decoupling of the (s)leptons which occurs 
when neglecting their Yukawa couplings, the full set of Boltzmann equations reduces to 
Eqs. (El), dHj), and ( ITOl) with the replacements Py6,yr5 ^mh ^ 0, Qi ^ 2((5 + T), coinciding 
with the set described in Ref. |21| and which we use in our numerical computation of the 
BAU. 



^ We thank C. Lee for this observation and T. Liu for highlighting this issue in an ear Her version of this 
work. 
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There are several uncertainties built into our computation of the baryon asymmetry. In 
addition to those arising from the vev-insertion approximation, theoretical uncertainties in 
several other parameters associated with the phase transition such as the bubble wall thick- 
ness, velocity, and variation of the Higgs vevs A/3 can introduce (9(100 GeV) uncertainties 
in the curves of constant baryon density, similarly to the case of higgsino/gaugino sources 
(see [lol] and references therein). For concreteness, we consider conservative values for the 
wall velocity, v^^ = 5/100, and thickness, = 5/Tc. Additionally, non-resonant sources 
such as those computed in Refs. [4] and related work yield results for the BAU that can dif- 
fer significantly from the values computed using the vev-insertion approximation, especially 
away from the resonance. To take into account the uncertainties associated with a precise 
calculation of the BAU, we show on selected plots curves corresponding to 10 x Yobs^ as a 
rough upper bound on the stop source scenario, as well as curves of 0.1 x Yobs as a more 
conservative estimate of the BAU. As we will see, our conclusions hold across this wide range 
of uncertainty. 



A. Parameter Space 

The baryon asymmetry produced by stop sources depend on the masses rn^^ which are 
temperature dependent. Since the thermal masses are constant at a given temperature, we 
can equivalently investigate the potential of stop sources to produce the observed BAU by 
varying the values of the LH and RH stop SUSY breaking soft masses, Mq^ jy^. We vary 
Mjj^ over the range -100^ GeV^ < M?^ < 2000^ GeV^ which includes the so-called "hght 
stop scenario" for M?^ < (and multi-TeV Mq^)^ a region of the parameter space where 
the light RH stop provides the strongly first order phase transition needed for successful 
baryogenesis. We stress that away from negative values of some other mechanism is 
needed to generate a strongly first order phase transition. Several such mechanisms have 



been proposed 3^, l38l that are decoupled from the spectrum required for EWB and from the 
physical processes of interest for the present discussion. Thus, in evaluating the potential 
of stops, sbottoms, and staus for EWB, we consider only the strength of the CP-violation 
in each case and assume a strongly first-order EWPT generated by one of these other 
mechanisms. For the LH stops, we vary Mg^ over the range 100 GeV < Mg^ < 4000 GeV. 
Within these mass ranges there are regions where the choice of soft mass leads to negative or 
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zero mass squared for the lightest stop at T = for various values of the triscalar coupling 
and /x; we indicate these regions (along with more stringent constraints on stop masses 
from direct searches discussed in Sec. IIIIBI) on all of our plots. For the stops, we focus on 
tan /3 = 10. One should note that larger tan /3 yields larger SM-like Higgs masses along with 
more stringent EDM constraints. 

There are several other parameters whose values need to be fixed in order to calculate 
the BAU. We choose values for these parameters conservatively, bearing in mind the various 
constraints from Higgs mass measurements, stop searches, and EDM search null results as 
well as theoretical considerations such as the avoidance of color and charge-breaking vacua. 
In computing the baryon asymmetry, we take ttia = 200 GeV. For larger m^, the baryon 
asymmetry is reduced due to the dependence of A/3 on m^, which scales as A/3 ^ V'^a 
{4]. The gaugino soft masses are taken to be real, with Mi = 80 GeV, M2 = 250 GeV 
to ensure a light neutralino the lightest supersymmetric particle (LSP) while other 

gauginos are rather heavy. For the scenario we consider here the resulting baryon asymmetry 
and Higgs mass constraints do not depend sensitively on Mi, M2. The gluino soft mass is 
largely decoupled from the phenomenology relevant here, and was set to M3 = 10 TeV. For 
the higgsino mass parameter /x (which we take to be real, so that 0^ arises only from the 
phase in ^4^), we choose /x = 200 GeV, 1000 GeV to illustrate the behavior of the baryon 
asymmetry and the various constraints in these cases. Small values of /x suppress the BAU 
(c.f. Eq. ffT4l)). while large values can make the zero-temperature physical stop masses 
squared negative by making the off-diagonal components of the mixing matrix large, as well 
as strengthen the various EDM constraints. Similarly, we vary the magnitude of the trilinear 
scalar coupling \ At\ = 100 GeV, 250 GeV, 1000 GeV; larger values of \At\ also result in larger 
exclusions from EDM constraints. We typically consider the case of maximal CP-violating 
phase (j)t = 7t/ 2 to show the maximal extent of the EWB-compatible parameter space. We 
rely on this phase to produce all of the baryon asymmetry, setting all other CP-violating 
phases 0^ = to isolate the contribution from the stop sources to the BAU. Otherwise, 
non-stop sources such as those arising from gaugino-higgsino-vev interactions will further 
contribute to the baryon asymmetry and there will be additional contributions to the EDM 
constraints. Finally, all other triscalar couplings are taken to be zero, and all other sfermions 
in our analysis are taken to be heavy, rrisf 10 TeV. As shown in Ref. |32|], this effectively 
decouples them from the network of transport equations, since superequilibrium and Yukawa 
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rates that can transfer charge density between SM particles and their superpartners vanish 
for any of the masses much larger than the temperature. As a result, the densities {/} 
appearing in the transport equations for these sfermions (e.g. Qi) correspond entirely to an 
SM particle charge density, kj = ki. 

Using this spectrum, we calculate the baryon asymmetry generated by stop scattering 
off of the bubble wall and outline regions of the stop mass parameter space suitable for 
successful EWB in Figs. [2H3l regions consistent with the observed value of are shaded. 
The contours shown correspond to maximal CP- violating phase sin^t, while for smaller 
phase the baryon asymmetry is suppressed as are the EDM constraints. Decreasing | sin 0^1 
does not open up any additional parameter space for EWB. Several important features of 
the sources are shown in Figs. EHS From Eqs. [HI [151 the CP- violating source is resonant for 
m^^ ^ m:[^. This manifests itself as a resonance for ^ M? + {M^.j^ — M^.j^) ^ in 
the parameter space as shown. 

Also, there is an increase in the venerated BAU for M? < 0. This feature arises far 

n 

from the resonance □ but is straightforward to understand from the form of the CP- violating 
stop sources. The quantities wr^l^ Sr^l entering into S~^^^ depend on the physical masses 
m~^^ = M'^^^^ + M^.j^j^ and for M?^ ~^t-r^ the physical mass mi^ 0. In this 
regime, the Boltzmann distributions in the integrand for Si begin to diverge for = 0, 
^^(fc) ^1/ (e'^'~^^ — l) which corresponds physically to the abundance of nearly massless 
squarks produced in the thermal bath. We emphasize that M^^ < does not result in a 
tachyonic stop in the unbroken phase as long as M?^ > — M|,.^ (the thermal masses are 
0(100 GeV) so this is not an issue in the parameter space we consider). On the other hand, 
the zero-temperature stop mass eigenstates can turn negative in some of the parameter space; 
the corresponding regions are of course ruled out by direct searches for stops, corresponding 
to the black shaded regions in Figs. [2H3l 

Finally, we find that there are regions for which the produced baryon asymmetry switches 
sign. This eflFect arises due to the competition between Q and T densities in the expression 
for n^, Eq. f fTOl ): since Q and T densities carry opposite sign, when Mg^ is small, kq/kg 
decreases and the T contribution can win out and drive > 0. In Figs. [2H3l regions with 
Yb > are shaded green, while regions for which < are shaded blue. Since at present 
^ We caution the reader that away from the resonance, there may also be non-resonant contributions to the 
sources [4] which we do not consider here. 
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tan(3=10, A =250 GeV, |i=1000 GeV 



tanP=10, A =100 GeV, |i=1000 GeV 




FIG. 2. Regions of the stop soft supersymmetry breaking mass parameter space consistent with the 
observed value of the baryon asymmetry resulting from stop sources for ji — 1000 GeV, \At\ — 250 
GeV (Left) and \At\ — 100 GeV (Right). Regions shaded blue (green) correspond to Yb > Yobs 
with Yb < (Yb > 0) for maximal CP-violating phase. The dotted blue contour on the left marks 
the region that would be consistent with stop-sourced EWB if the vev-insertion approximation had 
underestimated Yb by a factor of 10 (we omit this curve in subsequent plots). On the left we also 
show, by the darker shaded regions, the parameter space compatible with 10 x the observed BAU, 
i.e. the allowed regions if the vev-insertion approximation overestimated Yb by a factor of 10. 
Black shaded regions are excluded by stop mass direct searches; regions to the left of the thick red 
line are excluded by LEP Higgs mass bounds in both cases. Current constraints on the electron, 
neutron, and ^^^Hg EDMs are represented by the black dashed- dot, dashed, and dashed- double- dot 
lines, respectively, with regions to the left of each line ruled out by null results; the projected future 
reaches for d^, d^, and dug measurements are shown in magenta (where applicable). In both cases 
here, both the d^ and dn future sensitivities lie above the plane shown. For the \At\ — 250 GeV 
case, the Mercury EDM future sensitivity also lies above the plane shown. 



the phase (j)t is not experimentaUy constrained, either region can lead to the appropriate 
overaU sign for the baryon asymmetry through an appropriate choice of 0^. 
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tan(3=10, A =250 GeV, |Li=200 GeV 



tan(3=10, A =1000 GeV, |i=1000 GeV 




M~ [GeV] M~ [GeV] 



FIG. 3. Same as Fig. but for \At\ = 250 GeV, fi = 200 GeV (Left) and \At\ = 1000 GeV, 
ji — 1000 GeV (Right). For \At\ — 250 GeV, the > curve falls beneath the black shaded region 
and future measurements of the neuron EDM are expected to probe all parameter space shown. For 
\At\ — 1000 GeV, the expected reaches of d^, dn, and dng future measurements lie above the plane 
shown here. 



B. Stop and Higgs Mass Constraints 

Having calculated the BAU resulting from stop sources, one should ask how the stop 
mass parameter space consistent with successful EWB confronts various other phenomeno- 
logical constraints. We consider three types of constraints on our EWB scenario: stop mass 
constraints from collider searches, Higgs mass bounds, and electric dipole moment search 
null results. 

Zero-temperature stop masses have been constrained by direct searches for superpartners 
at LEP and the LHC (for particle spectra relevant here) to be m^^ > 96 GeV [401]. This 
lower bound arises from considering stop decays to the LSP, typically assumed to be the 
lightest neutralino. With our choice of gaugino masses, the lightest stop ii is heavier than 
X? in all of the parameter space so that this lower bound on m^^ is applicable. There are 
several specific cases in which the stop masses might be more tightly constrained, however 
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for generality we consider this lower bound for our scenario. We calculate the physical T = 



stop masses using FeynHiggs [4j| for the choices of parameters discussed above and indicate 
m^^ < 96 GeV on our plots by the black shaded region. Increasing \At \ leads to larger regions 
of parameter space for which the lightest stop falls below the lower bound. This is because 
the triscalar coupling appears in the off-diagonal entries in the stop mass matrix and large 
values reduce the value of the smaller eigenvalue. 



The mass of the SM Higgs has been constrained by LEP to be mho > 115.5 GeV |42L l43| . 
We use FeynHiggs to calculate the mass of the SM-like Higgs to two-loop order and indicate 
the LEP bound by a thick red line on our plots. In addition to the lower bound from 
LEP, recent preliminary results from both ATLAS and CMS experiments have indicated the 
possibility of a SM-like Higgs with rrihQ ^ 125 GeV |4J|. However, for light stops and small 
\At\^ the corrections to uiho arising from diagrams with stop loops typically needed to increase 
the mass of the SM-like Higgs in the MSSM are suppressed, and we find no parameter space 
consistent with rrihQ = 125 GeV. For larger \At\^ the stop loop corrections can be enhanced 
and the Higgs mass can be pushed up to rrihQ ^120 GeV (which we indicate on the plot 
corresponding to \At\ = 1000 GeV, // = 1000 GeV with a red dotted line), however we 
find that rrihQ = 125 GeV is difficult to obtain for our choices of parameters. We note that 
additional field content, such as the inclusion of a gauge singlet in the superpotential in 
e.g. the NMSSM, which may be required to provide a strongly first order phase transition 
in these scenarios, can result in large contributions to m/^Q, even at tree level. Thus, our 
Higgs mass contours should not be taken as strict exclusions, but as illustrating the tension 
encountered in the MSSM between light third generation scalars and a heavy SM-like Higgs. 



C. EDM Constraints 

CP-odd couplings in the MSSM will generally give rise to electric dipole moments (EDMs) 
of elementary fermions, nucleons, and neutral atoms. To date, no EDM has been experi- 
mentally observed, with the most stringent limits having been obtained on the EDM of the 
neutral Mercury atom js^ {\dHg\ < 3 x 10~^^e cm), electron (via the YbF molecule) j45| 
{\de\ < 1.05 X 10~^^e cm), and neutron [46] {\dn\ < 2.9 x 10"^^e cm). The non-observation 
of these EDMs places powerful constraints upon the strength of the CP-odd sources used 
in EWB (for a discussion of the constraints relevant to Higgsino-Bino-Wino driven MSSM 
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baryogenesis, see, e.g.^ Refs. jsl, Q|) • On-going efforts could improve the sensitivity of EDM 
searches by up to two orders of magnitude (for a review, see, e.g.^ Ref. [47]), suggesting the 
future possibihty of even more stringent constraints or the observation of an EDM with a 
magnitude consistent with the requirements of MSSM EWB. 

In order to analyze the impact of the present and prospective constraints, we use the pro- 
gram CPsuperH |48| to compute the relevant EDMs under different scenarios. In particular, 
when CP- violation is generated entirely by the phase cpt^ the largest contributions to the 
relevant EDMs are generated by two-loop graphs that give rise to the Weinberg three-gluon 
operator (rf^) as well as to "Barr-Zee" graph contributions to the elementary fermion EDM 
(rfj ) and quark chromo-EDMs (rf^). In addition, four- fermion interactions are generated at 
one-loop order, though the effects of these operators are typically suppressed. 

Before proceeding, we note that there exists considerable theoretical uncertainty in the 
computation of EDMs of strongly-interacting and many-body systems. In the case of di- 
agmagnetic atoms such as ^^^Hg, the dominant contribution to the EDM arises from the 
nuclear Schiff moment induced by CP-violating nucleon-nucleon interactions. In general, 
the most important contribution to the latter is a long-range effect arising from single pion 
exchange, wherein one of the pion-nucleon vertices (ttNN) is CP-odJf and the other CP- 
even. The CP-odd nNN interaction can be induced by the Weinberg three-gluon operator, 
CP-odd four-quark operator, and/or quark chromo-EDM operator, though in the MSSM 
the latter contribution typically dominates [49]. The computation of the atomic EDM, thus, 
encounters several sources of theoretical uncertainty: the calculation of the CP-violating 
ttA^A^ vertices from the underlying CP- violating interaction; the computation of the nuclear 
Schiff moment that generally requires a scheme for nuclear model-space truncation; and the 
corresponding atomic physics computation of the induced EDM. 

The computation of the neutron EDM is clearly less susceptible to theoretical uncertain- 
ties, as only those associated with hadronic effects enter. Nonetheless, these uncertainties 
can be substantial for both the neutron and atomic EDMs. For example, recent work by 

n 

the authors of Ref. [50] utilizing the QCD sum rule technique suggests that the sensitivity 
of the neutron EDM to quark EDM and chomo-EDMs may be a factor of five smaller than 
given by earlier work [49] that provides the basis for the CPsuperH code. In the case of the 
nuclear Schiff moment contributions to the ^^^Hg EDM, the code has yet to take into ac- 
^ Technically speaking, the interaction is odd under parity and time-reversal. 
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FIG. 4. Neutron EDMs for M^^ = 800 GeV, tan/3 = 10, = 1000 GeV and \At\ = 250 GeV. Red 
denotes negative values. Left: the three independent calculations of the neutron EDM. Right: EDM 
subcomponents using QCD sum rules. By far the largest contribution comes from the down-quark 
chromo-EDM d^, followed by the down-quark EDM . 



count state-of-the-art many body computations [51^ that imply substantial differences with 
the many-body calculations using a simplified, schematic nuclear interaction on which the 
code is based. Consequently, we caution that the precise numerical results associated with 
the neutron and ^^^Hg EDM constraints given below should be taken with a grain of salt 
(we comment more on the impact of this uncertainty on our results below). To provide 
an indication of the kind of theoretical uncertainty one might expect, we show in Fig. 4(a)| 



computations of the neutron EDM using diflFerent approaches as discussed in Ref. [48'] (QCD 
sum rules, the chiral quark model, and parton quark model), though we rely only on the 
QCD sum-rule technique in our analysis. The QCD sum rule computations tend to give the 
largest EDM, leading to the strongest constraints. 

With these caveats in mind, we observe that the strongest constraint for stop sources 
comes from the Mercury EDM, for which, in turn, the quark chromo-EDMs provide the 
most important contribution as shown in Fig. [5l We also include the constraint from the 
electron EDM. For the scenario of interest here, the largest contributions to the EDM 
come from top and stop loops in Barr-Zee graphs. Note that, like the CP- violating sources, 
all EDMs are roughly proportional to |//||74^|. Therefore, increasing one of these parameters 
in order to get a model with sufficiently strong baryogenesis also tends to produce a model 
that is ruled out by EDM searches. 
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FIG. 5. Left: a breakdown of the Mercury EDM, using the same parameters as in figure [H Almost 
the entire contribution comes from the down-quark chromo-EDM (multiphed by a constant factor). 
Right: a further breakdown of the down-quark chromo-EDM. 



D. Results 



The various constraints are plotted along with the curves of constant for different val- 
ues of /X, 1^4^ I in Figs. [2H3l For \At\ = 250 GeV, /x = 1000 GeV, one finds that direct search 
constraints rule out light, nearly-degenerate stop soft-breaking masses, with the bound from 
LEP on rrihQ excluding portions of the parameter space away from the resonance. Addi- 
tionally, null results from searches for the electron and neutron EDMs rule out nearly all 
of the EWB-compatible parameter space except for the tip of the resonance (the expected 
reach of future EDM experiments are also included in magenta). The strongest constraint 
is that arising from searches for the ^^^Hg EDM, which rules out all of the viable parameter 
space, even excluding regions in which the stops produce 1/10 of the observed BAU. This 
remarkable result is due to the stringent limit on the Mercury EDM coupled with the rel- 
atively large chromo-EDM contributions of diagrams involving stop loops to dng- We find 
a similar landscape for the | A| = 100 GeV, /x = 1000 GeV case in Fig. O here the BAU is 
reduced relative to the | A| = 250 GeV case and direct searches rule out less of the parameter 
space because of the reduced mixing. The Higgs mass constraints are stronger due to the 
smallness of \At\ and the LEP bound alone rules out all of the available parameter space 
for stop-sourced EWB. But once again the ^^^Hg EDM constraint is the most significant, 
excluding all parameter space in which the stops produce even 1/10 of the observed BAU. 
Moving on to the \At\ = 250 GeV, fi = 200 GeV scenario in Fig. [3l one might imagine 
that smaller values of /i, \At\ may reduce the impact of the EDM constraints enough to 
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open up some of the parameter space for stop-sourced EWB. However, although the EDM 
constraints are weakened, the baryon asymmetry is also reduced and once again the ^^^Hg 
constraints rule out all available parameter space. Finally, increasing both \At\ and /x to 
1000 GeV (see Fig. [3j) yields larger regions excluded by direct stop searches (due to the large 
mixing) and weaker Higgs constraints, allowing one to push rrihQ up to 120 GeV in parts 
of the parameter space. The EDM constraints in this case are more stringent, again ruling 
out all available parameter space for stop-sourced EWB. We note that in addition to the 
scenarios shown in Figs. [2}l3] we could have also chosen a small value of fi and large \At\^ 
however in this case the EDM constraints are again very stringent and all the stop-sourced 
EWB parameter space is excluded; we omit the corresponding figures for brevity (note that 
in this case one can obtain larger Higgs masses). We have also verified that decreasing sin^^ 
does not open up any more parameter space for stop sources. 

Additionally, varying tan/3 and/or uia does not aflFect our results. Since the stop-sourced 
baryon asymmetry and EDM predictions for light stop contributions scale with ^ 1/ sin/3, 
they are both rather insensitive to changes in tan /3 (we return to the large-tan /3 regime in 
the following section). Increasing tua suppresses both the overall baryon asymmetry and 
the expected EDMs, but the BAU varies as l/m\ whereas the EDMs only vary as l/m^, 
so a heavier CP-odd Higgs provides stronger exclusions. Smaller values of uia can enhance 
the BAU up to about a factor of 4 (for tua = 100 GeV), but we have checked that this does 
not overcome the strong EDM exclusions in any of the cases considered. 

As discussed above, the computation of dng involves significant theoretical uncertainty, 
which could impact the strength of the above conclusions. However, even if the theoretical 
prediction for dng were in fact an order of magnitude smaller than the values used in Figs. [21 
El it would still be just as constraining as the electron-EDM, which by itself rules out virtually 
all of the parameter space with \Yb\ > Yobs- We expect the vev-insertion approximation 
to ot'er-estimate the produced baryon asymmetry, and so it is unlikely that even this large 
correction would in fact open up any additional space for stop sources in a more careful 
treatment beyond the approximations used here. Similar considerations hold for the sbottom 
sources in Sec. [IV] as well. 

Summarizing, in considering the various scenarios depicted in Figs. [2}l3l we find no vi- 
able parameter space for MSSM stop- driven resonant EWB consistent with Higgs mass, stop 
mass, and EDM constraints. Even conservatively estimating the result of the various uncer- 
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tainties of the calculations as increasing Yb by an order of magnitude does not open up any 
viable parameter space for stop-sourced EWB. The large experimentally excluded regions 
are primarily a result of the stringent EDM constraints, and particularly that of ^^^Hg. We 
have also verified that even e.g. a factor of ten decrease in the Mercury EDM limits does not 
open up any additional parameter space for the stops. It is difficult to see how one might 
circumvent these constraints to produce the correct baryon asymmetry through a scenario 
relying primarily on stop sources. 



IV. THE LARGE tan/3 REGIME: SBOTTOM AND STAU SOURCES 

We now turn our attention to the other third-generation scalars as a possible source 
for the observed baryon asymmetry. Since the CP-violating sources arising from triscalar 
interactions for the sbottoms and staus are proportional to their Yukawa couplings, for these 
sources to contribute significantly to the BAU, one must consider larger values of tan/3. For 
large tan/3, the sbottom and stau Yukawa couplings, ^5, are enhanced as 

rrib rur . . 

Vb = Vr = 0^ 20 

V cos p V COS p 

where v ^ 175 GeV is the Higgs vev at T = 0. In what follows, we take tan/3 = 40 so 
that the strength of the sbottom and stau CP-violating sources are effectively comparable 
to that of the stops. 



A. Sbottoms 



To compute the CP-violating source for resonant sbottom scattering off of the EWPT 
bubble wall, we make use of the Lagrangian in Eq. ([T]). The sbottom interaction Lagrangian 
differs from that of the stops by the replacements {t, Af^ yt^ Vu^ ^^d} j^, Vb^ '^d, and 
the relevant CP- violating phase is (j)^ = AigdaAi^b'') . One can use these replacements in the 
non-equilibrium field theory derivation for S-'^^ in Ref. {21] to obtain S^^^ . The resulting 
CP- violating sbottom source is given by 

Ncy^ 
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Im 
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(21) 



where the Sl^r-, ojl^r terms are as in Eq. ( ITSll with Mi Mi and Ti corresponding to the 
thermal widths of the LH- and RH-sbottoms. Notice that the coupling of the sbottom to 
the down- type Higgs vev manifests itself as an overall relative sign between S^^^ and S^^^ . 
The source enters into the same set of QBEs, Eqs. (l7j)-(fT2l). and since there is now a source 
for one must include the density Qi in the network of equations. Since we take the 
sleptons to be heavy and the SM leptonic Yukawa interaction rates are small compared to 
the corresponding rates for the quarks, we neglect Ty^ in our calculation of the sbottom- 
sourced baryon asymmetry. As a result, the equations for the densities L and R decouple 
from the full set of QBEs; the relevant set of Boltzmann equations to solve is then given by 
Eqs. (CD, (ED, (Ej) and f fTOl) with the replacement Vy^ 0. In terms of the relevant charge 
densities, the left-handed fermionic charge density in Eq. (|6D simplifies to 

nL = ^Q + 2Qr. (22) 

We note that since Ir is heavy, the right-handed stops and tops are no longer in superequi- 
librium. This manifests itself in the contributions to the Yukawa and relaxation rates in- 
volving the RH stops vanishing, while the density T in Eq. ([7D corresponds entirely to a SM 
charge density. In the parameter space we consider the (s) bottoms are in superequillibrium 
everywhere except the kinematically disallowed region for the RH rate, and we proceed anal- 
ogously to the (s)top case by assuming ji^^ ^ = /x^^ ^ superequilibrium in all the parameter 
space when computing the baryon asymmetry. 

We calculate Yb following the spectrum outlined in Sec. IIII A[ only now with 100 GeV< 
< 2000, Mfj^ = 10 TeV, At = 0, and varying \Ai,\ = 100, 250, 1000 GeV. As for 
the stops, we assume that a strongly first order phase transition is generated from some 
mechanism other than the light stop scenario. The resulting sbottom-sourced BAU is plotted 
in Figs. [MZl where regions compatible with the observed asymmetry are shaded blue (we 
find no sign change in Yb for the sbottoms with our choices of parameters). The resonance 
in the CP- violating source is again apparent. 

In Figs. [MZl we also show the lower bound on the sbottom mass from direct searches, 
^ 89 GeV j4o| and contours of constant SM-like Higgs mass. The LEP bound on rrihQ 
is satisfied in all regions of parameter space considered. Since the mass of the RH stop is 
heavy, rrihQ receives larger contributions from stop loops compared to the stop-source case 
and one can easily push the Higgs mass up to ttz/^q ^ 120 GeV, however larger values are 
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FIG. 6. Regions of the sbottom soft supersymmetry breaking mass parameter space consistent with 
the observed value of the baryon asymmetry resulting from sbottom sources for ji — 1000 GeV, 
\A})\ — 250 GeV (Left) and \A})\ — 100 GeV (Right). Regions shaded blue correspond to Yb > Yobs 
for maximal CP-violating phase. The curve corresponding to an overestimate of Yb by a factor 
of 10 falls within the black shaded regions, which are excluded by sbottom mass direct searches. 
Red lines are iso-contours of the SM-like Higgs mass ruh^; the LEP bound is satisfied in all of the 
parameter space shown. The current constraint on the neutron and Mercury EDMs are represented 
by the black dashed and dashed- double dot lines, respectively, with regions to the left of each line 
ruled out by null results. For the \Ab\ — 250 GeV case, the current de bound falls beneath the 
shaded region while the current constraint on the Mercury EDM rules out all of the parameter 
space shown. Future EDM measurements of de, dn, and dng cltc expected to definitively probe well 
beyond the shown parameter space. For the \Ai)\ — 100 GeV case, the electron EDM current bound 
falls beneath the black shaded region, while future EDM measurements of de, d^, and dug will again 
probe all of the parameter space shown. 



more difficult to obtain with our choices of parameters (as with the stops, these should not 
be taken as strict exclusions). The EDM constraints for the sbottom sources are similar to 
those for the stop sources, but they receive a tan/3 enhancement, and thus the constraints 
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FIG. 7. Same as Fig. but for 1^1 = 250 GeV, fi = 200 GeV (Left) and \At\ = 1000 GeV, 
ji — 1000 GeV (Right). For both the \A})\ — 250 and \A})\ — 1000 GeV cases, the future reach of 
electron, neutron, and Mercury EDM measurements is expected to probe the entire parameter space 
shown. For the \Ai,\ = 1000 GeV case, the current Mercury EDM constraint already rules out all 
of the parameter space shown. 



tend to be more stringent. 

The behavior of the produced baryon asymmetry and the various constraints in Figs. [MZl 
is quahtatively similar to that for the stop-source case: increasing 1^451 or /x leads to larger 
regions compatible with the observed baryon asymmetry but strengthens the various EDM 
constraints. We note that since the sbottoms have down- type couplings to the Higgs, the 
roles of Ai, and /x in the mass matrix for the T = sbottoms are reversed relative to the roles 
of At and /x for the stops, and as a result, the exclusions from direct searches are primarily 
sensitive to /x for the large value of tan/3 chosen here. In addition to the cases shown in 
Figs.lMZl we have verified that scenarios for sbottom-sourced EWB with large 1^451 and small 
/X are also solidly ruled out by the current ^^^Hg EDM constraint. We have also checked 
that decreasing the strength of the CP-violating phase opens up no additional parameter 
space for the sbottom sources (it potentially could have, as the slope of EDM constraints 
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on the shown parameter space is different from that of B AU isolevel curves) . Consequently, 
taking shottom mass, Higgs mass, and EDM constraints into account, we find no regions 
of the shottom mass parameter space consistent with the observed value ofYs: as for stop 
sources, current EDM constraints imply that shottom sources alone cannot explain the BA U 
m the context of SUSY EWB. 



B. Staus 



FinaUy, we consider the case where the observed baryon asymmetry may have arisen 
primarily from CP- violation in the stau sector of the MSSM. For large values of tan 



can become enhanced as per Eq. f l20l ). From the Lagrangian, Eq. ([T]), and following Ref. j2l| . 
we can proceed in parallel to the calculation of Eq. ( |2T]) for the CP-violating stau source 
S^^^ with the replacements {6,745,^5} {f^A^^y^}^ yielding 
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and with the appropriate replacements in the definitions of Eq. f fTSl) for the LH and RH 
staus. The relevant CP- violating phase is now (j)^ = Arg(/xA^6*). The source Eq. f l23l ) enters 
the full set of QBEs, since for large tan /3 all third-generation Yukawa rates should be taken 
into account. The left-handed fermionic charge density is given by 

nL = Q + 2Qi + (24) 
kL 

where ki is the fermionic contribution to k^. Note that unlike for quarks, only the third 
generation LH density L contributes to since there is no generational mixing for leptons 
and we have neglected the first- and second- generation leptonic Yukawa couplings. We have 
verified that the staus and taus are in superequilibrium everywhere except in kinematically 
disallowed regions, so we proceeded as before, assuming /j^tlr = I-^tl r computing Yg. 

For the spectrum we again proceed in parallel to the analysis outlined in Sec. IIII Al 
with the appropriate replacements for the staus, again assuming a strongly first-order phase 
transition, either from the light stop scenario or some other mechanism (in calculating the 
BAU and constraints we assume a heavy RH stop). The resulting slepton-sourced baryon 
asymmetry is shown in Figs.[8H9]for various values of l^rL M maximal CP- violating phase 
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tanp=40, A^=250 GeV, |i=1000 GeV tanp=40, A^=100 GeV, |i=1000 GeV 




M~ [GeV] M~ [GeV] 

FIG. 8. ^5 in Figs. \MBi but for stau sources. For the \Ar\ = 250 GeV case, the dashed blue lines 
correspond to constant-Ys curves computed for a factor of ten larger thermal stau widths. In this 
case, the expected reach of future measurements will probe all of the parameter space shown. For 
the \Ar\ — 100 GeV case, neutron and Mercury EDM bounds fall beneath the black shaded region. 
In both cases the expected reach of future dng measurements is nearly degenerate with the current 
bound from measurements of d^ and is not shown. 



regions of the stau mass parameter space compatible with successful EWB are shaded 
blue. The resulting baryon asymmetry is strongly peaked near the resonance. This is because 
the thermal widths of the staus in the plasma, which enter into 8l,r in the denominator of 
S^^^ ^ are much smaller than those for the squarks, Tq^t — 0.5T, Vl^r 0.003T As 
a result, successful stau-sourced EWB requires nearly degenerate f^, tr] from Figs. [8H9]we 
find \Mi^ — M^^l ^ 100 GeV to produce the observed value of Yb for smcj)^ = 1- 

The results shown in Figs. [8H9] demonstrate that the resonance supplied by the small 
thermal widths of the staus present in the denominator of Eq. [23] can overcome the suppres- 
sion effect of the resonant relaxation rate F^^. This can be understood by noting that the 



overall baryon asymmetry scales with 



2l\ ^ Sf^^/\/T\^ so although both the source and 



relaxation rates are resonantly enhanced by the small widths, the asymmetry will tend to 
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tanp=40, A^=250 GeV, |i=200 GeV tanp=40, A^=1000 GeV, |i=1000 GeV 




M~ [GeV] M~ [GeV] 

FIG. 9. ^5 in Fig. O/or different values of A^, fi. For the \Ar\ = 1000 GeV case, future electron 
and neutron EDM experiments will probe all parameter space shown. In both cases the expected 
reach of future dug measurements is nearly degenerate with the current bound from measurements 
of de and is not shown here 



increase with decreasing widths. Also, the strong resonance in the denominator of S^^^ can 
overcome the Boltzmann suppression in the numerator for stau soft masses up to ^ 1 TeV in 
most cases. PhysicaUy, this corresponds to the very efficient production of chiral charge by 
a relatively small abundance of staus in the plasma. These results hinge on the small values 
Tlji = 0.003, which we take from Ref. and which were computed for tan/3 = 15. 
One might expect the widths to be enhanced for the larger values of tan /3 we consider here, 
since e.g. the otherwise negligible Yukawa decay f Hr can become important in this 
regime, yielding an enhancement of T^^ji from this extra decay channel of a factor of order 
2 at most. Also, T^^r are not necessarily equal, due to the differing hypercharges in the 
decays f tB^ as well as the SU{2) decay to rW which can be open for our choices of 
gaugino masses. To demonstrate the sensitivity of our results to the precise values of the 
thermal widths, we include on the LHS of Fig. [8] the curves of Yobs, 10 x Yobs calculated for 
an order-of-magnitude larger widths, F^.^^^ = 0.03T, which we expect to over-estimate the 
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uncertainty in Tl^r associated with these considerations. Even this factor of 10 increase in 
^L,R admits a significant amount of parameter space compatible with stau-sourced EWB. 
We thus expect that our overaU conclusions are rather insensitive to the details entering into 
a more precise determination of the stau thermal widths, however we encourage the Reader 
to keep the above caveats in mind when interpreting our results. 



In Figs. [HUH] we also plot constraints from direct searches for staus, rrif^ ^ 82 GeV |40|], 
which display the down-type dependence on and /i similar to that of the sbottoms. We 
also show iso-level contours of constant Higgs mass for the \Ar\ = 250 GeV, fi = 1000 GeV 
case. We omit these curves for the other plots since the exact values of the SM-like Higgs 
mass in each case are sensitive to the details of e.g. squark and gluino masses which do not 
impact the determination of Yb in slepton-sourced EWB. Two significant features emerge: 

(1) From Fig. [8] we see that one can achieve a Higgs mass in this scenario consistent with 
the hints from ATLAS and CMS, rrihQ ^ 125 GeV. In contrast to our analysis of the 
stop and sbottom sources wherein we found no viable regions of parameter space for 
rrihQ ^ 120 GeV, we are able to easily obtain a heavier Higgs mass for the stau source 
case since we are free to consider heavy squarks which contribute large loop corrections 
to rrihr, . Also note that the excess events observed in the H ^ channel with respect 

n n 

to general MSSM expectations |4J] could favor a scenario with light staus |35| . 



(2) Alternatively, one may also obtain the correct baryon asymmetry from stau sources 
along with a light right-handed stop; since there is no CP- violating phase in At^ the 
EDM constraints will not be affected, however the large loop contributions to rrihQ will 
be lost. 

We also consider EDM constraints on the stau-source scenario in Figs. [SHHl In order to 
generate chromo-EDMs, one needs a CP-violating phase that couples to (s)quarks. When the 
only phase is 0^, the chromo-EDMs disappear. Consequently, both the neutron and Mercury 
EDM constraints are much weaker in this scenario, while the electron EDM is the relevant 
one. The electron EDM in the case of CP violation in the stau sector entirely stems from a 
single Barr-Zee graph with a stau loop. From Figs. [SHHl we see that the lack of chromo-EDMs 
opens up large sections of parameter space, allowing for viable baryogenesis while satisfying 
the experimental constraints. Future EDM experiments are expected to probe all of the 
parameter space available for stau-mediated EWB. Note that since the primary constraint 
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on stau sources is the electron EDM, the available parameter space is rather insensitive to 
the theoretical uncertainty in the calculation of dng'- an order of magnitude ixnder-estimate 
of the mercury EDM would make its constraints comparable to that of the electron EDM. 
This picture holds even for smaller values of the CP-violating phase; we find that one can 
produce the correct BAU and still satisfy the various constraints for sin 0^ ^ 10~^ in most 
cases considered. 

Summarizing the results of Figs.lHHHl conclude that it is possible to produce the observed 
baryon asymmetry with CP violation in the stau sector only, for nearly degenerate staus and 
300 GeV ^ ^L3E3 ^ 1-2 TeV, depending on the magnitudes of the stau triscalar coupling 
and ji, while satisfying EDM and direct search constraints. This scenario can also naturally 
accommodate an SM-like Higgs mass rrihQ ^ 125 GeV for heavy squarks, or possibly a strongly 
first order phase transition via the light stop scenario for light mi^ . Should future searches 
for electron, neutron, and Mercury EDMs yield null results, all scalar sources in the MSSM 
will be ruled out for resonant EWB. 

V. DISCUSSION AND CONCLUSIONS 

The exploration of the electroweak scale is closing in on supersymmetric models in which 
the baryon asymmetry is produced at the electroweak phase transition. In this study, we 
have focused on the possibility that the CP violation necessary for successful EWB is found 
in the third generation sfermion sector, a scenario considered only for stops in the past. 
Here we have studied, for the first time, CP violation in the third generation down-type 
sfermion sector, i.e. sbottoms and staus, as a source relevant for baryogenesis. Also, we have 
quantitatively addressed the question of whether or not sfermionic CP violating sources are 
compatible with current constraints on the size of the electron, neutron and atomic electric 
dipole moments. 

The main findings of this study are 

(1) Neither the stop nor the sbottom sector are viable options to account for the bulk of the 
observed baryon asymmetry of the universe: two-loop Barr-Zee diagrams contribute to 
the chromo-EDM of the down quark to a level that is ruled out by current constraints 
of the Mercury EDM across the entirety of the parameter space where stop or sbottom 
sources could source a large enough amount of baryon asymmetry. Moreover, stop- 
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and/or sbottom-mediated EWBG is disfavored by indications of mho ^ 125 GeV, 
though present Higgs search constraints on the CP- violating sources are not nearly as 
decisive as those arising from EDMs. 

(2) The stau sector (where no chromo-EDMs are produced) has milder constraints from 
EDMs, and hence can be responsible for producing the net left-handed chiral charge 
density needed to produce, via weak sphaleron transitions, the observed baryon asym- 
metry in the universe. It is also possible in this case to achieve rrihQ ^ 125 GeV or a 
light RH stop as needed for a strong first order phase transition, but not both. Due to 
the relatively small stau thermal widths, however, this scenario of "slepton-mediated" 
electroweak baryogenesis requires almost degenerate staus, with masses between 300 
GeV and 1.2 TeV, depending on the size of the stau triscalar coupling and /x. This 
scenario also requires large values of tan/3. 

While, from the standpoint of requiring a strongly first order phase transition, electroweak 
baryogenesis in the MSSM is being conclusively tested by ongoing searches for a light stop 
and the Higgs jisl, the present results provide important, complementary information on 
the nature of the new sources of CP violation, the second key ingredient for successful 
supersymmetric electroweak baryogenesis. Future results from the LHC and from the next 
generation of EDM searches are therefore expected to yield an increasingly sharper, if not 
definitively clear, picture of whether or not the electroweak scale is related to the generation 
of the observed baryon asymmetry. 
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